Light beams with general direction and polarization: 
global description and geometric phase 



R Nityananda 1 and S Sridhar 2 

1 TIFR Centre for Interdisciplinary Sciences, 21, Brundavan colony, Narsinghi, 
Hyderabad 500 089, India; National Centre for Radio Astrophysics, TIFR, 
Pune 411 007, India 

2 Raman Research Institute, Sadashivanagar, Bangalore 560 080, India 
E-mail: rajaram@ncra.tif r . res . in, ssridharOrri . res . in 

Abstract. We construct the manifold describing the family of plane monochromatic 
light waves with all directions, polarizations, phases and intensities. A smooth 
description of polarization, valid over the entire sphere S 2 of directions, is given through 
the construction of an orthogonal basis pair of complex polarization vectors for each 
direction; any light beam is then uniquely and smoothly specified by giving its direction 
and two complex amplitudes. This implies that the space of all light beams is the six 
dimensional manifold S 2 x C 2 , the Cartesian product of a sphere and a two dimensional 
complex vector space. A Hopf map (i.e mapping the two complex amplitudes to the 
Stokes parameters) then leads to the four dimensional manifold S 2 x S 2 which describes 
beams with all directions and polarization states. This product of two spheres can be 
viewed as an ordered pair of two points on a single sphere, in contrast to earlier work 
in which the same system was represented using Majorana's mapping of the states of a 
spin one quantum system to an unordered pair of points on a sphere. This is a different 
manifold, CP 2 , two dimensional complex projective space, which does not faithfully 
represent the full space of all directions and polarizations. Following the now-standard 
framework, we exhibit the fibre bundle whose total space is the set of all light beams 
of non-zero intensity, and base space S 2 x S 2 . We give the U(l) connection which 
determines the geometric phase as the line integral of a one-form along a closed curve 
in the total space. Bases are classified as globally smooth, global but singular, and 
local, with the last type of basis being defined only when the curve traversed by the 
system is given. Existing as well as new formulae for the geometric phase are presented 
in this overall framework. 



PACS numbers: 42.25. Ja, 03.65.Vf, 02.40.-k, 42.81. Gs 



Submitted to: J. Phys. A: Math. Gen. 



Global description of light beams and their geometric phase 



2 



1. Introduction and summary 

Bertolotti pQ formulated the evolution of linear polarization as light traverses a space 
curve in an inhomogeneous but locally isotropic medium in the geometrical optics limit. 
He concluded that the electric vector is parallel-transported with respect to a connection 
derived from a conformally flat metric, where Euclidean distances are scaled by the local 
value of the refractive index. Rytov [2] derived this evolution law independently from 
a WKB treatment of Maxwell's equations, and also expressed it as a phase difference 
per unit length between the two circular polarizations, proportional to the torsion of 
the space curve. Vladimirskii |3J brought out the following geometrical implication: 
polarization vectors live in the tangent plane to the sphere of directions and undergo 
parallel displacement as the direction changes. This implies that after the tangent vector 
to the curve returns to its original value (e.g. after one turn of a helix), the polarization 
rotates by an angle equal to the solid angle enclosed by the closed trajectory of the 
tangent vector on this sphere. 

Pancharatnam [I] , in the context of novel interference patterns shown by absorbing 
biaxial crystals, formulated the phase which now bears his name, equal to one-half of the 
solid angle traversed on the Poincare sphere which represents polarization states. The 
work of Berry j5] on the phase change of a quantum state, evolving adiabatically as the 
Hamiltonian describes a closed path in a parameter space, and its later generalizations, 
provides the natural framework in which to discuss this class of optical situations; see 
e.g. [HI [7] for reviews of early work. The Berry or geometric phase depends on the path 
traversed in the parameter space (i.e the sphere of directions or of polarizations) but 
not on the rate of traversal. 

Bhandari [8j [9] , Hannay pi)] and Tavrov et al [H] , among others, treated geometric 
phases under the simultaneous evolution of direction and polarization. This can be 
viewed as occurring in a four dimensional space, whose global structure is naturally of 
interest. Bhandari and Hannay worked with the four dimensional ray space (the space 
of physical states) of a spin one quantum system. Such a ray space is constructed from a 
three dimensional complex vector space (Hilbert space) by identifying vectors differing 
only in normalization and phase; it is denoted by CP 2 , two dimensional complex 
projective space. The spin one description is indeed a faithful mapping of the elliptical 
orbits traversed by the electric fields of the light beams being considered, which are the 
same as the orbits of a three dimensional isotropic harmonic oscillator. However, two 
problems prevent CP 2 from being the four dimensional space that faithfully represents 
all directions and polarizations. One is that a given ellipse traversed by the electric 
vector in a plane in (real) three dimensional space can correspond to two opposite 
directions of propagation. This could be resolved by doubling CP 2 , using one copy for 
each sense of circular or elliptic polarization. However, another problem arises on the 
boundary between these two regions: an electric field which is linearly polarized along 
a given direction can belong, not just to two but, to an entire circle of directions of 
propagation in a plane perpendicular to it; i.e. the subspace of linear polarizations in 
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the spin one model is two dimensional, whereas it should really be three dimensional. 
Thus the correct global nature of the space is worth investigating. 

In this paper we first determine the six dimensional manifold which faithfully 
represents light beams with all directions, polarizations, phases and intensities. As a by- 
product we give a direct demonstration that the four dimensional space that represents 
all directions and polarizations is S 2 xS 2 , the Cartesian product of two spheres, which is 
distinct from CP 2 . As discussed below this implies more than the known fact that each 
of direction and polarization forms a sphere, when considered separately. In principle, 
the total space of interest could have been a non-trivial, twisted S 2 (polarization) fibre 
bundle over S 2 (directions) as basejj] Given a direction of propagation with unit vector k, 
two real orthogonal unit vectors can be chosen as a pair of linearly polarized basis states 
in the tangent plane. Then the complex electric vector, of any general polarization, can 
be specified by affixing complex amplitudes to the basis vectors. These two complex 
amplitudes contain information about polarization, phase and intensity. Hence the six 
dimensional space of all light beams can be thought of as a fibre bundle with base 
S 2 the sphere of directions, and fibre C 2 a two dimensional complex vector space. A 
smooth extension of the locally constructed real, linearly polarized, basis vectors to all 
directions k G S 2 runs into the obstruction that it is impossible to construct a smooth 
non- vanishing vector field on S 2 . We show below that this obstruction can be overcome 
by choosing a complex basis. 

In section 2 we exhibit such an orthogonal pair of complex basis vectors for each 
direction which is globally smooth on S 2 . Section 3 brings out some implications of 
this construction. The existence of such a globally smooth basis implies that the fibre 
bundle is trivial. Therefore the six dimensional manifold is S 2 x C 2 . It then follows from 
the standard transition from complex amplitudes to the Stokes parameters ( the Hopf 
map) that the four dimensional manifold which describes beams with all directions and 
polarizations is S 2 x S 2 . Given a general state of polarization with a general direction, 
one can pick a unique point k on the sphere of directions and, using this basis, assign a 
unique point s on the Poincare sphere. 

The geometrical phase, discovered by Berry in the context of adiabatic evolution 
of a quantum-mechanical system, was restated in the language of modern differential 
geometry by Simon |12| . as the (an)holonomy of a U(l) connection on a fibre bundle. 
Samuel and Bhandari [13] offered a general setting for the geometric phase, allowing 
for comparison of the phase between any two nonorthogonal states in Hilbert space. In 
section 4 we follow the formulation of [13] and construct a fibre bundle, whose total 
space is the set of all light beams of non-zero intensity S 2 x C 2 \ {0} , with base space 
S 2 x S 2 . We use Pancharatnam's method of comparing phases of different polarizations 
and Vladimirskii's notion of parallel displacement to derive a U(l) connection on the 
fibre bundle, and give a basis-independent definition of the geometric phase. In section 5 
we use the global basis discussed above, and write the geometric phase as the sum of two 

| Incidentally, CP 2 cannot be represented as such a bundle. 
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phases; one which resembles the Pancharatnam phase |4j and the other which resembles 
the phase of Rytov and Vladimirskii [3j. Each of these terms depends on the global 
basis used; but their sum, the geometric phase, is independent of basis. The result of 
Hannay [TO] follows simply and concisely from one particular choice of global (but not 
smooth) basis. The bases used in [21 El El EH EH EE] are all local, not global. We define 
general local bases and discuss two particular cases of interest, in section 6. 



2. Construction of a globally smooth basis 

We begin with the pair (0 , <p) of unit vectors along the orthogonal directions of 
increasing [9 ,(p) . This basis is ill-defined at the north and south poles but smooth 
everywhere else. The following sequence of basis pairs are constructed and their 
properties listed below, with the final one being the required globally smooth pair. 

I m — cos V 9 ~ sxntpip , 

In2 = simpO + cos (ftp. (1) 

This linearly polarized pair is smooth at the north pole since the transformation 
compensates for the counterclockwise rotation of the (8 , (p) directions there. However, 
this doubles the clockwise rotation of (6 , ip) at the south pole, and hence this basis is 
singular at 9 = n . The linearly polarized pair 

I si = + cos (p 6 + sin<^<£, 

ls2 = —sm.ip6 + cos (2) 

is, similarly, smooth at the south pole and singular at the north. For each of the three 
linear bases {(0 , cp) ; (l N1 , l N2 ) ; (I si , ^2)} there is a corresponding circular basis. These 
are defined by 

cr = -j= {0 + iip) , 

cl = (0 - itp) ; (3) 



c NR = —7=(lm + Urn) = exp[+i(p]c R 



V2 
1 



c N l = -7=(lm ~ iM = exp[-iy3]c L ; (4) 



c S r = + = exp[-iip]c R , 

csl = (I si - iisa) = exp [+i<p}c L . (5) 



Global description of light beams and their geometric phase 



5 



Notice that all the circular bases have well defined state of polarisation everywhere but 
have phase singularities. The (c R , c L ) basis has singular phase at both the poles, while 
the N and S sets are smooth at N and S respectively, but are singular at the opposite 
poles, S and N. 

The globally smooth basis is given in terms of (c R , cl) by 

9i = Aexp[+i V ]{cos 2 (9/2)c R - sin 2 (9/2) c L } , 

g 2 = Aexp [-i<p] {sin 2 (6/2) c R + cos 2 (9/2) c L ) , (6) 

with the normalization factor A = [cos 4 (0/2) + sin (9/2)^ . The (g 1 ,g 2 ) basis is 
constructed to tend to (c^ R , c^l) at the north pole and (— c$l , csr) at the south pole 
(note the switch in order in the latter case), and is hence smooth everywhere on the S 2 
of directions. Expressing (g ± , g 2 ) in terms of (0 , ip) 

g 1 = B exp [+i<p] {cos 9 9 + 'up} , 

g 2 = Bexp[-i<p]{0 - icos9cp} , (7) 

with the normalization factor given by B = [1 + cos 2 9] . This form brings out 
a simple physical interpretation - it is readily seen that the basis vector g l is right 
circularly polarized at the north pole of the sphere of directions. As we move southwards, 
it becomes elliptically polarized with axial ratio cos 9 and the major axis along the <p 
direction, thus linearly polarized on and along the equator and reversing the sense of 
traversal of the ellipse in the southern hemisphere and becoming left circular, with a 
unique well-defined phase that is independent of <p at the south pole, this is exactly 
the pattern of polarization and phase (though not intensity) which would be radiated 
by an electric dipole rotating in the positive sense in the xy plane, viewed as a function 
of direction. The second basis vector g 2 is left circular at the north pole and elliptical 
in the northern hemisphere with its major axis along the direction; it reverses sense, 
retaining the same major axis, as it evolves to become right circular at the south pole 
- one can think of it as the electric field radiated by a magnetic dipole rotating in the 
negative sense in the xy plane. 

Complementing the above physical realisation, we offer another simple way of 
appreciating the existence of globally smooth nonvanishing complex vector fields on 
a sphere. Both sin 9 6) and sin 9 ip are smooth real vector fields on the sphere, vanishing 
at the south and north poles and winding once around each to make up the total index 
(Euler characteristic) of two. Such fields can of course be constructed for any pair of 
antipodal poles by rotation. One then manufactures a nonvanishing complex vector 
field simply by ensuring that these antipodal nulls are not the same for the real and 
imaginary part. It is straightforward to normalize one such field, and to construct the 
other member of theth orthonormal pair by interchanging real and imaginary parts 
and changing the sign of any one of them (this corresponds to turning the major axis 
of a general elliptical polarization by 7r/2 and reversing the sense, steps which preserve 
smoothness). Our choice, (J7J), has the real part of g l vanishing on the y axis (of direction 
space) and the imaginary part vanishing on the x axis. The above choice of basis is 
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far from unique. As discussed in the next section, any other global basis is related to 
(g 1 , g 2 ) through a U (2) valued field on the direction sphere. So all smooth choices of 
bases belong to a single homotopy class. 



3. Consequences of the global basis 

The globally smooth basis provides us with two complex vector fields gi(k) and g 2 (k) 
on the base S 2 of directions. The complex electric vector of any light wave propagating 
along the k direction can be written as 

E = z x g x {k) + z 2 g 2 (k), (8) 

where Z\ and z 2 are any two complex numbers. Hence the six dimensional space of light 
beams with all directions, polarizations, phases and intensities is S 2 x C 2 , which is the 
(Cartesian) product of a sphere and a two dimensional complex vector spacejjjl This 
means that once the smooth basis field is given, the triplet (k , Z\ , z 2 ) , giving direction 
and two complex amplitudes, specifies the plane light wave completely. The intensity 
of the beam is 

E*-E = \ Zl \ 2 + \z 2 \ 2 , (9) 

because <?*• g l — g\- g 2 — 1 and g\ - g 2 = g 2 - g x = . The same beam can be described 
in a rotated basis which is also globally smooth. We write the new basis as 

(g[(k) , g' 2 {k)) = , fl 2 (*0) M- X (fc) , (10) 

where M(A;) is an U(2) valued field on S 2 . The new complex amplitudes are given by 
corresponding U(2) rotations, (z[,z' 2 ) T = M(k)(zi , z 2 ) T , which preserves the form of 
the intensity given in (Q. 

Beams of non-zero intensity have E*- E > . The space describing these beams is 
T = 5 2 xC 2 \ {0} , where <C 2 \ {0} is the set of all (zi, z 2 ) except for the origin (0, 0) . For 
a given direction k, two points, [z\ ,z 2 ) and (z[ ,z 2 ), have the same polarization state 
(or ellipse) if z 2 /z[ = z 2 j z\ . Hence the natural complex coordinate on the Poincare 
sphere is ( = z 2 jz\ , which includes the point at infinity when z\ — and the phase of 
z 2 is irrelevant to the polarization state. This relationship is defined by the Hopf map 
from C 2 \ {0} to S 2 , given by 

/:C 2 \{0} -> S 2 , {z 1 ,z 2 )^z 2 /z 1 . (11) 

Polar coordinates, (9 p ,(p p ), can be defined by ( = tan(6 , p /2) exp(iy9 p ) . Then the 
Stokes parameters (in units of the intensity) for fully polarized light correspond to 
the components of the unit vector s : 

z\z 2 + Z 2 Z\ 



z\zi + z%z 2 



sin 9 P cos if p ; 



§ The fact that making the tangent bundle on the sphere complex also makes it trivial, is well-known 
to some mathematicians we have consulted. 
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z\z 2 — Z\Z\ 
Z\Z\ + Z 2 Z 2 



) 



sin 6 P sin (p p ; 



Z\Z\ — Z 2 Z 2 



cos 6 P . 



(12) 



Z*Zi + Z 2 Z2 



The ordered pair of unit vectors (k ,s) determines uniquely a light beam with direction k 
and polarization state labelled by s . Therefore the four dimensional manifold combining 
the direction and state of polarization of light is £ = S 2 x S 2 the product of two 
spheres. It should be noted that a given point s may correspond to different polarizations 
states for different k. If we want to know the parameters (axis ratio, handedness and 
orientation) of the polarization ellipse, we require the information contained in the basis 



4. Parallel— transport of E 

4-1. Fibre bundle describing light beams of non-zero intensity 

Here we construct a fibre bundle with total space T and base £, with projection given 
by using the Hopf map ffTTj) . The projection 7r g requires a choice of basis, which could 
be any of the global bases which are related to each other through the U(2) valued field 
M(fc), as discussed earlier; this is indicated by the subscript g. Then there are global 
coordinates (k,zi,z 2 ) on T. Under the action of 7r g , there are global coordinates on 
(fc, C = z-i jz\) on C . We define parallel-transport on T, using Pancharatnam's principle 
[I] of comparing phases of light beams of different polarizations through interference. 
Since we are dealing with two beams in infinitesimally close directions, we first have to 
carry out parallel trannsport on the sphere of directions to bring the second beam to 
the same direction as the first, consistent with Vladimirskii [3] The result may be stated 
simply: E is said to be parallel-transported when its local change in phase is zero; i.e. 
the projection of each complex field onto the previous one in the sequence is in phase 
with it. The notion of parallel-transport in T is formalized as an U(l) connection on 
the fibre bundle. In the construction and use of this connection we follow the framework 
presented by Samuel and Bhandari [13]. In the case treated by them, the total space 
is the set of (nonvanishing) states of a Hilbert space, which is not so here, but the 
construction goes through without any change. 

The fibre bundle (T , vr g , C , T , A) consists of the following elements: 

(i) The total space is the set of all nonzero intensity light beams T = S 2 x C 2 \ {0} , 



(ii) The base space C = S 2 x S 2 is the set of all light beams of all directions and 
polarizations. 

(iii) The (typical) fibre is J 7 ~ C \ {0} , the complex plane minus the origin. This set of 
non-zero complex numbers is contractible to the unit circle S 1 G C, so that J 7 is the 
set of all complex numbers to which a unique phase (modulo 27r) can be assigned. 



field ( gi (k),g 2 (k)). 



defined by 



T = {(k,E) I k-E = 0, E*-E>0}. 



(13) 
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(iv) The projection 7r g from the total space to the base depends on the global basis 
chosen; the subscript g refers to the global basis (flfi(fe) , 9 2 (k)) °f ©• Any of the 
other global bases discussed earlier can also be used. For instance, 7r g / can refer to 
the projection map when the basis (g[(k) , g 2 (k)) is used. Henceforth we assume 
that a choice of the global basis has been made. Then E is given by the ordered 
pair of complex numbers (z±, z 2 ) as in (jSj), so that (k, z lt z 2 ) are global coordinates 
on T . The Hopf map of ffTTl) can be used to define the projection we want: 

7r g :T->£, (k,z 1 ,z 2 ) H> (k,( = z 2 /z x ) . (14) 

(v) The inverse image, 

n~\k,0 = {(k,z l ,z 2 )eT\z 2 /z l = C}, (15) 

is the fibre over the base point (k, Q , which is diffeomorphic to the (typical) fibre 
J 7 ~ C \ {0} . Moving along any fibre corresponds to changing both the intensity 
and phase of a given polarization state. Changes of phase are of interest to us, so we 
imagine following the changes on the unit circle to which the fibre is contractible. 
Therefore, at fixed direction and intensity, the 3-sphere of polarizations and phases 
is realized as the total space of a non-trivial Hopf bundle with phase as S 1 fibre 
over the Poincare sphere S 2 of polarizations (see Nityananda [17] for an elementary 
treatment of this). 

(vi) Given any curve (k(£) ,E(£)) in the total space T, with a real parameter £, the 
tangent vector to the curve is (dk/d£ , dE/d£ = U{£)). We define 

A, = ^> (16, 

as the contraction of a connection one-form on the tangent vector, and will refer 
to it henceforth as the connection for brevity. Physically, Ai measures the rate of 
change of phase along the curve, and the sign is so chosen that a positive value 
corresponds to phase lagging as £ increases. We note some salient properties of An 
below: 

(a) An is a linear functional on the tangent space at any point of T . The definition 
given above is independent of basis. 

(b) Of the two components, dk/d£ and U(£), only the latter contributes to Ag . 
This is because Ai measures the rate of change of the phase of E with respect 
to the parameter £ , as is evident from ( |22l) ; phase changes are insensitive to 
the principal-curvature, dk/d£, of the space curve corresponding to the curve 
k(£) in direction space. 

(c) Gauge transformations move E along the fibre: E{£) — > p(£) exp [ia(£)]E(£) 
with p(£) > . Since Ai measures rate of change of phase with respect to £ , 
it is independent of p(£) which measures change in magnitude. Hence Ai is a 
U(l) connection on the fibre bundle, which transforms inhomogenously as: 

Ae ~* Ae+ d£- (17) 
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4-2. The U(l) connection 

Interference between two fully polarized beams E\ and E 2 traveling in the same 
direction gives, 

|| E 1 + E 2 || 2 = E* 1 -E 1 + E* 2 -E 2 + 2Re{E{- E 2 } . (18) 

Pancharatnam jl] noted that E\ and E 2 , referring to different non-orthogonal 
polarization states, can be said to be 'in phase' when, for fixed intensities of the 
individual beams, the total intensity of the superposition is maximized. Thus the phase 
difference is taken to be zero when {E\- E 2 } is real and positive. More generally, he 
defined the phase lag of E\ with respect to E 2 to be the argument of E\ - E 2 . We (and 
the earlier authors cited) apply the same concept to the case where the electric field 
vectors lie in infinitesimally close planes rather than in the same plane. The component 
of the electric field vector associated with the new direction perpendicular to the original 
plane is first order in the angle between the planes, but it does not contribute to the 
scalar product. 

Let the curve, (k(£) , E(£)), in the total space T represent a light beam that is 
taken around the curve, k(£) , in direction space. The electric field E(£) is defined to be 
parallel-transported when it changes least. By this we mean that E(£) is in phase 
with E(£ + d£) , while k(£)- E(£) = k(£ + d£)- E(£ + d£) = 0. Therefore parallel 
transport implies that lm{E*-U} = 0. We derive the U(l) connection and discuss 
its consequences, using methods given in Samuel and Bhandari |13j: 

(i) A curve in T on which Ai = is defined to be a horizontal curve. A horizontal 
curve in T is said to be a horizontal lift of its projection onto the base. 

(ii) The horizontal lift of a closed curve in £ may be open in T; i.e the initial and final 
points lie on the same fibre but have different phases and magnitudes. The phase 
difference is the geometric phase, formally described as the (an)holonomy of the 
U(l) connection on the base. 

(iii) Let {k(£) , E(£)) be a horizontal curve in T such that, as £ varies from to 1, the 
light beam returns to the same direction and polarization state. This means 

fc(l) = fc(0), E{1) = XE(0), (19) 

where A is a non zero complex number. The initial and final states are on the fibre 
over the same base point, and the final lags the initial by an amount Arg{A}, which 
is the geometric phase. 

(iv) Define a closed curve Ch by completing the circuit in T by joining (k(l) , E(l)) to 
(k(0) , E(0)) by a segment running along the fibre. It can be verified that 

7 = I A e d£ = -Arg{A} + 2mr , (20) 

where n is the winding number of the (vertical) segment running along the fibre. 
This is because, along the horizontal curve = , and all the contribution comes 
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from the vertical segment which is equal to — Arg{A} + 2mr jjj] Note that 7 as 
defined in fl20l) is minus the conventional geometric phase, which is the phase of 
E(l) with respect to E(0) . Henceforth we refer to 7 as the geometric phase. 

(v) Ch projects to a closed curve C in the base £ = S 2 x S 2 , which is the locus of 
points (k(£) , ((£)) , where 

fc(i) = Mo), C(i) = C(o). (21) 

The integral in ( 1201) is along a specific closed curve Ch in the total space T which 
projects to a closed curve C on the base. However, it is gauge invariant, because 
Ai transforms according to ( fTTl) . Therefore 7 is a property of C (up to an integral 
multiple of 2ir), making the geometric phase factor exp [ry] a property of the closed 
curve C. 



5. The geometric phase 

Let C be a closed curve on C = S 2 x S 2 with parameter i that varies from to 1 . 
Let C be a closed path in T, obtained by lifting C . Gauge invariance implies that the 
geometric phase is given by, 

A lit =P^l M . (22) 

Thus Ai measures rate of change of phase, along the curve C , with respect to the 
parameter i . This expression also makes it clear that the geometric phase is independent 
of the basis used. Just like in (]2"U]1. the value of the integral in (J22"]) can differ by integral 
multiples of 2tt , when it is evaluated along paths in T that are both lifts of C but not 
homotopic to each other. Once a basis is chosen, C acquires coordinates: these could 
be the ordered pair (fe, £), or (k, s) where s is the Stokes vector. 




5.1. Globally smooth basis 

We write the electric field in the global basis (<7i(fe) ,g 2 (k)) as 

E{£) = z x {l) 9l {k{&)) + z 2 (£) g 2 (k(£)) , (23) 
and define the state vector, 

* - { 3 ) • 

which is an element of Af ~ C 2 \ {0}, the nonzero states of a two dimensional Hilbert 
space H . When z\ and z 2 return to their original values along the closed path C, their 
arguments ^1 and ^2 can change only by integral multiples of 2tt : 

6(1) = 6(0) + 27H7T, 6(1) = 6(0) + 2n 2 7r, (25) 

|| Since the fibre C \ {0} is contractible to S , curves from a given initial to a final point on the same 
fibre can differ by an integral number of windings around this S . 
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where n\ and n 2 are integers. Since z 2 / z 1 = tan(6 l p /2) exp(iy? p ) , we have <p p = (£ 2 _ £i)) 
and the change in (p p along the path C is: 

^ p (l) = ^p(O) + 2(n 2 -m)7r, (26) 

so (n 2 — rii) measures the winding number of the azimuth (p v on the Poincare sphere, 
and is hence a property of the curve C on the base C . For a given path on the base 
- i.e keeping (n 2 — rii) fixed — n\ counts winding on the fibre, and different values 
correspond to homotopically distinct lifts C of the path on the base. We find it useful 
to keep the explicit dependence on both these integers (which is commonly left implicit) 
in the discussion which follows. 

Substituting ( 123|) in (TIB"]) , we obtain, 

_ Im{M(d/dl)M} (07 , 

where if is a Hermitian operator — henceforth referred to as a 'Hamiltonian' — which 
defined on H (and hence A/"), whose matrix elements are given by, 

1 (_, dfl, da* 



"« = iA 9 "'-iir-ir-e>)- (28) 

Of the two terms on the right side of (l2"7j) . the first can be written (upto a total derivative) 
in terms of coordinates on the Poincare sphere. Direct computation gives 

ht y^ = i(i-«.^ + d{„ (29) 

The second term of (|27p can be thought of as the expectation value of H in the state 
l^) . It depends on both the rate of change of the basis along the curve, as well as the 
polarization state. 

The total geometric phase can now be written as 

7 = 7p + 7rv, (30) 
the sum of two quantities. The first is a Pancharatnam-like (hence P) phase: 

7p = % — m — de 



\ j> (1- cos6» p )dv? p + f d^i 
Jc Jc 



= |fi p + 2n l7 r, (31) 

where fl p is the solid angle on the Poincare sphere swept by the shortest meridional arc 
that connects the north pole to the Stokes vector s , and the contribution of 2niir comes 
from the total change in £i over the closed path cj^j] Even though the value of 7p is 

% Note that J7 P changes discontinuously by ±4-7r when the contour crosses the south pole, contributing 
an additional amount ±^4-7r = ±2n to 7p, which is harmless because it does not change the geometric 
phase factor exp [ry]. 
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basis-dependent, its form is invariant, and equal to |O p + 2ni7r . We refer to the second 
term, 

as a Rytov-Vladimirskii (hence RV) phase, since they derived it for a situation when 
there are changes of direction but polarization remains circular. The phase 7rv is more 
general than [21 [3] in that it also keeps track of changes of polarization. The elements 
of the Hamiltonian H can be computed by substituting (JJJ) in ( 128]) : 

sin 2 

H u d£ = -H 22 di = — -dip, 

1 + cos 

H 12 di = H* 21 dt = - ^'^^ [^n 26 d^ + 2i d6] . (33) 

Note that the use of the global basis gives a Hamiltonian which is smooth everywhere: 
terms singular at the poles are multiplied by factors which vanish there. Using this is 
( |32l) . it is straightforward to calculate the RV phase as a line integral over the closed 
curve C on the base: 

/ sin 2 9 [cos 9 P — cos 9 sin 9 V cos (o? p — 2<p)] 

7RV = f c TT^d ^ 

+ /^gsing p sin(y -2y) d6) _ 

Jq 1 + COS 2 9 

The total geometric phase, 7 = 7p + 7rv, is obtained by adding the expressions in (I3T]) 
and (El]). However, the separation should not be taken literally, since both 7 P and 7rv 
depend on the basis chosen. It is only their sum 7 which is independent of basis. 

We can also take account of geometric phases due to lumped elements introduced 
in the path of the light ray, causing abrupt changes of polarization. Let the state vector 
change by a finite amount, from to when £ changes by an infinitesimal amount 
di . Since the change in direction is smooth, the change in k is also infinitesimal; the 
integrand of fl32|) being finite, the contribution to 7rv during this jump in polarization 
is of order d£, and can be neglected. However, 7p gets a finite contribution, which is 
given by the geodesic rule of [13]: the phase change is given by adding a line segment 
AC to the integral of (!3T|) . such that AC* is a geodesic curve in M, connecting and 

. This projects to a geodesic, AC , which is a great circle on the Poincare sphere, 
connecting the initial and final Stokes vectors, s and s' . Hence, phase changes due to 
multiple lumped elements encountered by a light ray can be accounted for, by adding 
geodesic arcs on the Poincare sphere (along which evolution is continuous) to the curve 
C in C. 



5.2. Global but singular basis 



Whereas the expression in (134"]) for 7rv is smooth in the global basis, it is difficult 
to interpret physically, because the global basis is variable on the direction sphere. 
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However, we can now utilize the freedom to use a non-global basis set for computation 
of any quantity whose existence is guaranteed by a globally smooth basis. Using 
the circular basis (cnr , c^l) of (HI) means dealing with the familiar azimuthal phase 
singularity at the south pole of our coordinate system. Let the state vector in this basis 
be \CN) = (z[ , z' 2 ) T ■ It is straightforward to calculate the Hamiltonian Hcn as, 

Hcn = (l-cos0)^f * _° \ , (35) 

Using this in (1271) . the connection one-form is, 

lm{(CN\(d/d£)\CN)} (CN\H CN \CN) 
1 "' (CN\CN) (CN\CN) 

= \ (1 — cos 9 P ) dip p + d£i + cos# p (l — cos 9) dip , (36) 

which is, as expected, singular at the south pole (6 = 7r) of the direction sphere. 
Integrating, the geometric phase is given by the compact expression, 

7 = <f> A e d£ = |fip + 2n X T\ + & cos P (1 - cos 9) dip , (37) 
Jc Jc 

as the sum of a P phase and an RV phase. We have seen that the former has an invariant 

form. The RV phase has the following structure: the basis being circular, cos 9 p = V, 

the Stokes V parameter which is a measure of the degree of circular polarization; and 

(1 — cos 9) dip is equal to the change in the solid angle traced by the shortest meridional 

arc connecting the north pole to the contour on the direction sphere (the CN basis being 

singular at the south pole, we assume that the contour does not traverse the south pole). 

Note that the physical meaning of the separation into P and RV phases, is specific to the 

chosen (i.e. CN) basis. We will see more examples of such basis-dependent separation 

in the next section. 

The work of Hannay |10| . using the Majorana representation for spin one, is 
formulated in terms of vectors in R 3 (ordinary three space), and gives a concise one- 
term formula for the geometric phase (equation (16) in [IQ]). This formula is manifestly 
intrinsic, i.e it is expressed in a form which makes it clear that it is independent of 
the basis chosen in polarization and direction spaces. We now place this result in 
our framework. Working in the CN circular basis as earlier, we write the complex 
electric field as E = ziC NR + z 2 c NL . We now make the gauge choice £i = —<p p /2; since 
<P> P = (£2 — £1), this implies that £2 = V 9 p/2, so that 

z l = \E\ cos(0 p /2) exp[-i^ p /2] , z 2 = \E\ sin(fl p /2) exp[+iy? p /2] . (38) 

It can be verified that the zero of phase is at the tip of the major axis of the ellipse: i.e. 
at time t = the real electric field vector points along the major axis which itself makes 
an angle of ip p /2 with l N i . This choice of phase, apart from being undefined at both the 
poles of the Poincare sphere, runs into a discontinuity of tt in phase as we go around the 
sphere once in the azimuthal direction, since the ellipse turns by 180° and the zero of 
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phase is transferred to the opposite end of the major axis. Therefore it is necessary to 
introduce a branch cut which could be chosen along a line of constant longitude joining 
the two poles. However, neither of these affect the application of this gauge in the small. 
The formula for the connection can be obtained by setting £ x = — <^ p /2 in (|36|) : 

Aidi = — icos^pdyjp + cos P (1 — cos9)dip = —cos9 p d/3, (39) 

where 5(3 = [^d<p p — (1 — cos#)d<^] has the physical interpretation of an intrinsic or 
covariant change in the azimuth of the major axis, with the second term allowing for 
the rotation of the (Ini,In2) frame under parallel displacement along the ip direction, 
because there is no rotation when we move along 6 . Equation (139]) is equivalent to 
equation (16) of [TO], so we see that the results of that paper follow from a particular basis 
choice in our framework, and are not dependent on using the Majorana representation. 
The cautionary remarks following equation (16) in [TU] follow from the existence of the 
branch-cut in this basis (which the curve C may traverse). 



6. The geometric phase in local bases 

The bases used in J2J H2 El dU [EH H5] are all of a local character; i.e. they are 
specified once the curve in parameter space is given, but not in advance. The bases 
also differ from each other, so that the expressions for 7p and 7rv, obtained by different 
authors, take different forms. It is straightforward to write the geometric phase when 
a local basis, (ei(£) , e 2 (£)) , is used. The local basis is defined on the curve, k(£) by, 
(ei(£) , e 2 {£)) = (g^k^)) , g 2 (k(£))) R~ l (£) , where R(£) is a £7(2) valued field which is 
defined on the curve in direction spaceQ This basis is orthonormal: e*- e.\ = eJj- e 2 = 1 
and e*- e 2 = e 2 - e± = . The new complex amplitudes are (wi , u> 2 ) T = R(fe)(zi , 2; 2 ) T . 
Then the complex electric field can be written as, 

E(£) = wi(£)ei(£) + w 2 (£)e 2 (£). (40) 

Using this in ( 1271) . the connection one-form is: 

= — W) — + ~mT' ^ (41) 

where {w\{£) ,w 2 (£)) T = \4>{£)) = R(£)\ip(£)) and K is the new Hamiltonian K, whose 
components are, 

Kl > 2i V 1 d£ d£ V • 1 j 

6.1. Local circular basis 

A local circular basis satisfies the usual orthonormality conditions, c\ - c\ = c 2 - c 2 = 1 
and c\ • c 2 = clj • C\ = . The Hamiltonian 

dci 
'~d£ ~ 



K ■ — — 
2i 



dc 2 


I 




~dl _ 







(43) 

+ The local basis can be extended, if necessary, to a tube of trajectories surrounding the curve k(£) . 
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is diagonal because, under a parallel displacement the circular basis vectors do not 
mix; each can only acquire a phase0 It is straightforward to calculate the connection 
one-form: 

V 

A t dl = i(l-cosfl p )d^ p + d& + — (cj-dci - c* 2 - dc 2 ) , (44) 

where V = cos8 p is the Stokes V parameter, measuring the degree of circular 
polarization along the curve. We can express the second term of (1441) in terms of 
properties intrinsic to the space curve. Begin with the linearly polarized basis used by 
Rytov [2]: (n(l), 6(0) , where n(Z) is the normal and b(l) is the binormal to the space 
curve described by the light ray; the parameter I measures length along the space curve. 
Define the Rytov circular basis as, 

J R {1) = ^ WO + 16(01 , 4(0 = ^ MO - , (45) 

and use the Frenet-Serret formulae to get, (2i) _1 (c\*- dc\ — c"- dc r 2 ) = —r(l)dl, where 
t(0 is the torsion of the space curve. Then the geometric phase is, 

7 = j)_A t dl = |fi p + 2nm - j>V{l)r{l)dl, (46) 

as before the sum of two phases. The second is the line integral of the product of Stokes 
V and torsion r , a property intrinsic to the space curve described by the light ray. For 
circularly polarized light, when V = ±1 , we encounter the integral of the torsion over 
a closed curve which equals (2tt — O), consistent with the known result for this case. 
The expression given above in (146p resembles ( 1371) for the geometric phase in the CN 
basis (which is global but singular). However, the former has the advantage of using 
quantities intrinsic to the curve. Circular bases are used in [14j [15l [16] . 

We address the concern that the normal and binormal are undefined for a straight 
section, or at a point of inflection, of our space curve. If the curve happens to be entirely 
straight, any orthogonal pair of real vectors in the transverse plane would suffice. The 
case of interest is therefore when part of the space curve has a well-defined curvature 
and normal. Working with the corresponding path on the sphere of directions, we see 
that each point represents a unit tangent k to the space curve. The unit tangent to the 
path on the sphere is just the normal n to the original space curve. A straight section 
is then just a single point, but the curve must continue in order to return to its starting 
point and, for a smooth continuation, there is no difficulty in specifying n, and hence 
the binormal 6 = fcxn, in a continuous manner. 

6.2. Parallel-transported basis 

Any basis, chosen at an initial point 1 = 0, can be evolved by parallel displacement. A 
useful property of an initial circular basis is that it always remains circular, acquiring 
only phase as it is evolved along the curve by parallel displacement. We take the initial 
basis at I = to agree with the Rytov circular basis, but undo the phase change which 

* We have already encountered an example of this in ()35|) . for the circular basis (c^ , cl) ■ 
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it undergoes by parallel displacement. Therefore, from (l4"6"j) . we see that the pair of 
circular states defined by, 

Pl (l) = exp [+i /„' r(J') dl>] c\(l) , p 2 (l) = exp [-i jj r(J') <U'] c r 2 (Z) , (47) 

is a parallel-transported circular basis; [8j [11] use bases like these. Unlike the bases 
used so far, the initial (at I = 0) and final (at I = 1) basis vectors are not the same. 
The vectors acquire a phase of minus (for right circular) and plus (for left circular) the 
solid angle ft in direction space: 

Pl (l) = exp [-iQ] Pl (0) , pa(l) = exp [+ift] p 2 (0) . (48) 

When the complex electric field, 

E(l) = w 1 (l)p 1 (l) + w 2 (l)p 2 (l), (49) 
returns to the same fibre, E(l) = XE(0) , the initial and final amplitudes are given by 

u>i(l) = Aexp[+ift]tz>i(0) , w 2 {l) = A exp [-ift] w 2 {0) . (50) 
As discussed in item (v) of section 4.2, we have the freedom to evaluate the integral on 



a convenient closed curve in the total space, so we can set A = 1 in (1501) . If 6 and £2 be 
the arguments of W\ and w 2 , we have 

6(1) = 6(0) + ft + 2 ni n, 6(1) = 6(0) - ft + 2n 2 vr, (51) 

where n\ and n2 are integers. Then, 

^ = exp[-i2 Si ]^), (52) 

implies that the curve C op on the Poincare sphere (associated with the complex 
coefficients (101,102) in the local basis) is, in general, not closed; its end points have 
equal latitudes but different longitudes: 

P (1) = P (O), ^p(l) = <^ p (0) - 2ft + 2(n 2 - ni )n. (53) 

We emphasize that this auxiliary Poincare sphere, associated with the local basis, should 
not be confused with the sphere of polarization in a global basis, which occurs in the 
definition of the S 2 x S 2 base space; on this latter sphere the path is of course closed. 
When ( fj^U is substituted fflUfy , the Hamiltonian K vanishes; this is the characteristic 
property of the parallel-transported basis. Hence the connection one-form is given by 
just the Pancharatnam term, 

dl = hKyiw} , = , (1 _ cos e) d + d6 (54) 

(010) 

Integrating, and using f loTj) . the geometric phase is: 

7 = [Adl = \ I (l-cos0 p )d^ p + [6(1) - 6(0)] 

J C J C p 

= iftjf + ft + 2n l7 r, (55) 
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where is the solid angle on the Poincare sphere, swept by the shortest meridional 
arc joining the north pole to the Stokes vector, as the latter traverses the open curve 
C op whose end points are given by ( 1531) . Equation ( 1551) is a compact and pleasing form: 
apart from the usual 2ni7T term, the geometric phase is given as the sum of one-half of 
the swept-out solid angle on the Poincare sphere and the solid angle on the direction 
sphere. For right circular light, the solid angle on the Poincare sphere = 0, and the 
geometric phase is simply (Q + 2n 1 7r), as it should be. However, the use of this formula 
calls for a limiting process when applied to the case when the polarization remains left 
circular throughout, since there are multiple geodesies connecting the north pole to the 
south pole. The degeneracy is lifted by introducing an infinitesimal ellipticity. From 
( 153]) we see that the initial and final longitudes on the Poincare sphere are well-defined 
and differ by [—2ft + 2(n 2 — ni)7r], and the solid angle swept out by the meridian 
is twice this amount. Then ( 1551) gives the correct phase, equal to (—ft + 2n 2 n) for left 
circular polarization. 

7. Conclusions 

The manifold describing the family of plane monochromatic light waves with all 
directions, polarizations, phases and intensities is seen to be S 2 x C 2 , the (untwisted) 
Cartesian product of a sphere and a two dimensional complex vector space. The four 
dimensional manifold, which describes beams with all directions and polarizations, is 
S 2 x S 2 , the Cartesian product of two spheres. The fibre bundle, whose total space is the 
set of all light beams of non-zero intensity S 2 xC 2 \ {0} , with base space S 2 x S 2 , has 
a (7(1) connection which allows a basis-independent definition of the geometric phase. 
When expressed in a general global basis, the geometric phase decomposes into the sum 
of two terms, one of which resembles the Pancharatnam phase, and the other the phase 
of Rytov and Vladimirskii. Each term depends on the global basis used, but their sum, 
the geometric phase, is independent of basis. 

The existence of a globally smooth description allows the use of non-global bases 
for computations. The geometric phase assumes a compact and simple form in terms 
of the polar coordinates on the Poincare and directions spheres. The final result for the 
geometric phase in Hannay's work [10J follows, in our framework, from a particular choice 
of basis (singular at the poles and along a cut joining them). Two examples are given of 
oppositely circularly polarized bases, which are not defined globally, but only over the 
path traversed by the system. The first follows Rytov in using the normal and binormal 
to the space curve, and the second uses parallel displacement, both giving interesting 
forms for the geometric phase. It is intriguing that one seems to be forced to choose 
between global smoothness and physical transparency: the smoothest form with which 
we begin is the most opaque, whereas the parallel-displaced local basis (which does 
not even return to the same value at the end of a closed path in direction/polarization 
space) gives a compact and pleasing result. 

We conclude with the global conceptual difference between our framework and the 
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widely used spin one representation. The latter admits an inner product between any 
two states, not just two nearby ones: the total space is C 3 , the space of normalized 
vectors is S 5 , and the space of states is CP 2 . We find that the space of normalized 
vectors for this problem is S 2 x S 3 , which is not a vector space and does not have a 
globally defined inner product. An example bringing out this global distinction is a pair 
of oppositely circularly polarized waves of unit intensity travelling in opposite directions: 
they have identical electric fields and hence would have a vector space inner product 
with a modulus of unity, even though they are not the same state by any means. The 
U(l) connection given above in (fl6|) — which determines the geometric phase accrued 
when a closed path is traversed in S 2 x S 2 — relies only on the inner product between 
two nearby electric fields. In mathematical terms, the set of all possible directions and 
their associated complex transverse electric fields is a vector bundle with a connection, 
which is not a vector space. We have presented a framework that describes the global 
structure of the space of directions and polarizations, and this allows us to view earlier 
work on the geometric phase for light in a unified way. 
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